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ABSTRACT 

An analysis of Lyman series lines arising from hydrogen-like oxygen and neon 
ions in the coronae of the active RS CVn-type binaries II Peg and IM Peg, 
observed using the Chandra High Resolution Transmission Grating Spectrograph, 
shows significant decrements in the Lya/Ly/? ratios as compared with theoretical 
predictions and with the same ratios observed in similar active binaries. We 
interpret these decrements in terms of resonance scattering of line photons out 
of the line-of-sight; these observations present the first strong evidence for this 
effect in active stellar coronae. The net line photon loss implies a non-uniform 
and asymmetric surface distribution of emitting structures on these stars. Escape 
probability arguments, together with the observed line ratios and estimates of the 
emitting plasma density, imply typical line-of-sight sizes of the coronal structures 
that dominate the X-ray emission of 10 10 cm at temperatures of 3 x 10 6 K and 
of 10 8 cm at 10 7 K. These sizes are an order of magnitude larger than predicted 
by simple quasi-static coronal loops models, but are still very small compared to 
the several 10 11 cm radii of the underlying stars. 

Subject headings: Radiative transfer — X-rays: stars — stars: coronae — stars: 
late type 



1. Introduction 

The Sun is the only star for which we can presently image coronal X-ray emission, 
and in order to understand how coronae on other stars might be structured we must resort 
to indirect means. One potentially powerful diagnostic of the characteristic size of X-ray 
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emitting regions is line quenching through resonance scattering. The escape probability of 
a photon emitted by a resonance line in a low density homogeneous plasma is dependent on 
the line-of-sight path length through the plasma region. Provided the plasma density and 
the abundance fraction of the ion in question are known, in the regime where the plasma is 
only marginally thick the comparison of strong and weak lines subject to different scattering 
losses can yield an estimate of the typical emitting region path length. 

Resonance scattering in the solar corona has been studied predominantly in the light of 
the strong (gf = 2.66) 2p 6 1 Sq — >■ 2p 5 3d 1 P± resonance line of Fe XVII at 15.02 A as compared 
to nearby weaker Fe XVII lines, though with controversial results concerning whether optical 
depth effects were seen or not (Phillips et al. 1996; Phillips et al. 1997; Schmelz et al. 1997; 
Saba et al. 1999). Studies of the same transition seen in different stars observed at high 
resolution (A/AA ~ 1000) by Chandra have also recently been presented by Phillips et al. 
(2001) and Ness et al. (2003). Both stellar studies failed to find evidence for significant 
optical depth. 

There are two problems with using the prominent Fe XVII soft X-ray complex for optical 
depth studies. Recently, it has been shown by Doron & Behar (2002) and Gu (2003) that 
the indirect processes of radiative recombination, dielectronic recombination, and resonance 
excitation involving the neighbouring charge states are important for understanding the 
relative strengths of Fe XVII-XX lines. Secondly, the coronae of active stars have been 
found to be Fe-poor by factors of up to 10 compared with a solar or local cosmic composition 
(e.g. reviews by Drake et al. 2003, Audard et al. 2003), reducing the sensitivity of Fe lines as 
optical depth indicators. Indeed, the spectral lines likely to exhibit the largest optical depths 
in the coronae of active stars are the Lyman a lines of hydrogen-like O and Ne — elements 
that are often seen to be enhanced relative to Fe (e.g. Drake et al. 2001, Audard et al. 2003). 

In this Letter, we present an analysis of the Lya to Ly(3 line strength ratios in the 
active binaries II Peg and IM Peg, and show that the Lya lines of Ne and O are significantly 
quenched. We use inferred optical depths to make the first direct estimates of the dimensions 
of coronal structures in active stars. 



2. Observations and Analysis 

Observations and relevant characteristics of the program stars are summarised in Ta- 
ble 1. All spectra were obtained by the Chandra High Energy Transmission Grating (HETG) 
in conjunction with the Advanced CCD Imaging System spectroscopic detector (ACIS-S), 
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and were downloaded from the Chandra Data Archive 1 . Spectra from multiple observa- 
tions (IM Peg and AR Lac) were combined, and positive and negative orders were summed, 
keeping HEG and MEG spectra separate. The same observations of HR 1099, II Peg and 
AR lac are described in more detail by Drake et al. (2001), Huenemoerder et al. (2001) and 
Huenemoerder et al. (2003), respectively. 

Resonance scattering of Lya and Ly/3 photons can be diagnosed by comparison of the 
measured Lya/LyS ratio with respect to the theoretical ratio. Spectral line fluxes were 
measured using the FITLINES utility in the PINTofALE 2 IDL 3 software package (Kashyap 
& Drake 2000). In the case of O VIII, the Ly[3 transition is blended with an Fe XVIII line 
(2s 2 2p 5 2 P 3/2 - 2s 2 2p 4 ( 3 P)3s 2 P 3/2 , A = 16.004 A). We estimated its intensity by scaling the 
observed intensity of the neighbouring Fe XVIII 16.071 A (2s 2 2p 5 2 P 3/2 -2s 2 2p 4 ( 3 P)3s 4 P 5/2 ) 
transition, that shares the same upper level, by the ratio of their theoretical line strengths 
(0.76) as predicted by the APED database (Smith et al. 2001). Uncertainties involved in 
this deblending procedure are negligible for II Peg (whose Fe XVIII is weak) and similar to 
or less than Poisson errors for the other stars for, e.g., an error of ~ 20% in the theoretical 
Fe XVIII ratio. He-like resonance lines of Ne IX and O VII were also measured in order 
to derive a temperature estimate from the Lya/He-like r ratio. Measured line fluxes and 
statistical errors are listed in Table 2. 

Observed and theoretical O and Ne Lja/Lj(3 ratios are compared in Figure 1. The 
observed ratios are shown at the temperatures at which observed H-like to He-like 2p-ls line 
intensities matched their theoretical (APED) values (see Table 3). While it is an approxima- 
tion to assume that the lines are formed at a single temperature, this is not critical because 
the temperature dependences of these Lya/Ly/3 ratios are not steep for T > 2 MK. Both O 
and Ne H-like lines are formed at temperatures significantly lower than those at which the 
theoretical Lya/Ly/3 ratios reach their asymptotic limits. 

The observed Lya/Ly/3 ratios are lower than the theoretical values for both the Ne X and 
the O VIII lines in the case of IM Peg, and for O VIII in II Peg. Huenemoerder et al. (2001) 
noticed similar discrepancies in II Peg but excluded optical depth effects on the grounds that 
this deviation was less than ~ 2a from the theoretical value; it appears, however, that this 
assessment was based on the high temperature asymptotic ratio of ~ 6.25 (in photon units). 
Our measured ratio for II Peg is instead more than 3<r lower than the expected ratio at the 
temperature of formation of the O VIII lines. In the case of IM Peg, the observed O VIII 
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ratio lies ~ 1.8a below the theoretical value, while the Ne X ratio in both HEG and MEG, 
are > 3a lower. For comparison with IM Peg and II Peg, we also present similar results 
for AR Lac and HR 1099: these ratios instead do not show significant departure from the 
expected optically thin values. 

Though uncertainties in the theoretical ratio are not included, these are not expected to 
exceed 10% based on good agreement with recent laboratory experiments (G. Brown, private 
communication; Beiersdorfer 2003). If we assume a 10% error in theoretical ratios, the Ne 
Lya/Ly/5 of IM Peg still departs by 3a for MEG (~ 1.7a for HEG), while the O ratios lie 
2.5a for II Peg and 1.3a for IM Peg below the theoretical curve. In Figure 1, we also plot the 
O VIII Lya/Ly/3 ratios obtained from the line fluxes measured by Raassen et al. (2002) from 
Chandra/LETGS and XMM- Newton/RGSl spectra of Procyon. In the Procyon spectrum 
Fe XVIII is not detected and these data provide further validation of the APED O VIII 
Lja/LjP ratio. 

These quantitative results are reinforced by a visual comparison of spectra in Figure 2, 
in which the Lya lines in II Peg and IM Peg are visibly weaker relative to Ly{3 that those of 
AR Lac and HR 1099. There do not appear to be plausible explanations for the discrepant 
ratios other than by the quenching of Lya relative to Ly/3 through resonance scattering 
within the emitting plasma. Intervening photoelectric absorption could cause similar effects 
but H column densities of order a few 10 21 cm 2 would be required — two or three orders of 
magnitude higher than inferred for these stars (Mewe et al. 1997; Mitrou et al. 1997). 



2.1. Path Length Estimate 

Our measured line ratios allow us to derive an effective optical depth r, and a typical 
photon path length within the emitting plasma. We use the escape probability, p(r), derived 
by Kastner & Kastner (1990), which for r < 50 can be approximated by (e.g. Kaastra & 
Mewe 1995) 

v(t) ~ . (1) 

n ; 1 + 0.43 r 1 ; 

This is the escape probability for line photons emitted at optical depths between and r, 

averaged over a Gaussian line profile due to thermal Doppler broadening, and it assumes that 

each scattered photon is completely lost from the line of sight. This probability is significantly 

larger than the e~ T transmittance of the simple absorption case, since it assumes emission 

over the whole line of sight through the plasma. The line center optical depth, r, can be 
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written (e.g., Acton 1978): 



-14 ™i A n H w \M 



r = 1.16 • 10- 14 • ^A z -^Xf\-n e £ (2) 
n ei n e V T 

for ion fraction n\/n e \ (from Mazzotta et al. 1998), element abundance Az, oscillator strength 
/, temperature T, electron density n c , atomic weight M, where nu/n e ~ 0.85, and £ is 
the total path length along the line of sight through the emitting plasma. In the above, 
neglect of non-thermal broadening should not be important: non-thermal velocities of, e.g., 
50 km s _1 — similar to the thermal velocities for O and Ne at their characteristic temperatures 
of formation — would result in £ being underestimated by a factor of 1.4. 

Electron densities were adopted from the survey of Testa et al. (2004; hereafter Paper I): 
for t(0 VIII) we assumed the n e derived from the O VII He-like triplet; the Ne IX triplet 
is strongly affected by blends with Fe lines and so for r(Ne X) we assumed the n e derived 
from the Mg XI He-like triplet that forms at a similar temperature. In the case of IM Peg, 
the low signal in the O VII intercombination and forbidden lines precludes a firm density 
estimate; we therefore assumed a value of 2 x 10 10 cm -3 , which is typical of values obtained 
for all the measurable spectra of RS CVns in Paper I. The assumed n e are listed in Table 3, 
while the adopted O and Ne abundances are listed in Table 4. 

In order to derive an estimate of the path length, £, we treated the fine structure 
components of Lya (1: 2 P 3/2 -* 2 S 1/2 at A = 18.9671 A and 2: 2 P 1/2 -+ 2 S l/2 at A = 
18.9726 A) separately, since their splitting (AA/A ~ 0.0003) is larger than the thermal 
width (~ a//cbT/M /c ~ 0.00013). The observed and theoretical intensities, J b s and J t h, are 
then related as follows: 



' Lya obs 



— h obs + h obs — 



th 



+ 



'2 th 



(3) 



1 + 0.43ri 1 + 0.43r 2 

Since for hydrogenic ions f 2 = fi/2 = 0.2776/2 (e.g. Morton 2003), I 2 th = h th/2, and 
71 ~ C(£) ■ fi (see Eq. 2), we obtain: 



i- Lya obs 
I Lya th 



+ • (4) 

1 + 0.43C(£)/i 1 + 0.43C(£)/i/2j v ; 

For Ly/3 the splitting of the components is smaller than the thermal width and the escape 
probability for LyS is given directly by Eq.l. The combined equation for C(£) was then 
solved to obtain the path length. 

The path length estimates, £ T) derived from the measurements are listed in Table 4. For 
comparison, we also list the stellar radii and loop lengths expected for a standard hydrostatic 
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loop model (e.g. Rosner et al. 1978, RTV hereafter), Lrtv, corresponding to the observed 
temperatures and densities. 

3. Discussion and Conclusions 

The discrepant O VIII and Ne X Lya/Ly/3 ratios found here for II Peg and IM Peg 
represent the first clear evidence of resonant scattering in coronal X-ray emission lines. The 
photon path lengths inferred from the observed ratios (Table 4) are: (1) about two orders of 
magnitude different from each other, reflecting the differences in the plasma densities found 
in Paper I for the characteristic temperatures of formation of the O VIII and Ne X lines; (2) 
very small with respect to the stellar radius; (3) an order of magnitude larger than the loop 
lengths derived from RTV model scaling laws. In the ratio £ t /Lrtv the n e terms cancel, so 
that this conclusion is completely independent of plasma density measurements. 

We note that the path lengths obtained from this type of analysis should be treated as 
lower limits because addition of Lyct photons by any optically thin emission or any scatter- 
ing into the line-of-sight are not taken into account. In the case of a uniform spherically- 
symmetric arrangement of emitting structures over the stellar surface, scattering into and 
out of the line-of-sight would be expected to compensate one another. The detection of 
resonance scattering implicitly suggests a non-uniform coronal distribution. 

That we find different optical depths in emitting regions on what are ostensibly similar 
active RS CVn-type binaries is puzzling. One possibility is that this is transient behaviour 
that might occur on all similarly active stars for a preferential arrangement of one or more 
coronal structures. We also note that, whereas both AR Lac and HR 1099 comprise two 
late-type stars of more similar type, II Peg and IM Peg have unseen companions of unknown 
spectral type, though how this arrangement might benefit larger coronal photon path lengths 
is not obvious. 

We thank the referee, Jeffrey L. Linsky, for helpful comments that enabled us to im- 
prove the manuscript. PT was partially supported by Chandra grants GO1-20006X and 
GO1-2012X under the SAO Predoctoral Fellowship program. JJD was supported by NASA 
contract NAS8-39073 to the Chandra X-ray Center, EED was supported by NASA grant 
NAG5-10872. GP and PT were partially supported by MIUR and by ASI. 



Table 1. Properties of Observed Stars and related HETG observations. 



Source 


Spectr.Type 


d (i) 


R*/R© 


P k ( 2 ) 
A orb 


Prot « 


L x [HEGf 


L x [MEG] a 


V b 

fx 


Obs ID 


^cxp 






[PC] 




[ksec] 


[ksec] 


[erg/sec] 


[erg/sec] 


[10 5 erg/cm 2 /sec] 




[ksec] 


II Peg 


K2V/.. 


42 


3.4/.. (») 


580 


580 


1.56-10 31 


1.76-10 31 


250 


1451 


42.7 


IM Peg 


K2III-II/.. 


96.8 


13/.. W 


2100 


2100 


2.75-10 31 


2.79-10 31 


27.1 


2527 


24.6 


IM Peg 












2.17-10 31 


2.30-10 31 


22.3 


2528 


24.8 


IM Peg 












1.86-10 31 


1.97-10 31 


19.2 


2529 


24.8 


AR Lac 


G2IV/K0IV 


42 


1.8/3.1 P) 


170 


170 


5.21-10 30 


5.60-10 30 


284 


6 


32.1 


AR Lac 












5.61-10 30 


6.30-10 30 


319 


9 


32.2 


HR 1099 


G5IV/K1IV 


29.0 


1.3/3.9 P) 


250 


250 


7.85-10 30 


1.05-10 31 


1020 


62538 


94.7 



a HEG range: 1.5-15 A; MEG range: 2-24 A. 

b X-ray surface flux from L x obtained from MEG spectra ^ 

I 



Note. — References: W SIMBAD database; ( 2 ) Strassmeier et al. 1993; ( 3 ) Berdyugina et al. 1998; ^ Berdyugina et al. 1999 



Table 2. Line flux measurements and Lya/Ly/5 ratios, with la errors. 



Source grating flux (10 6 photons cm 2 sec x ) Photon Ratio 

Nc X Ne IX O VIII O VII Fc XVIII Lya/Ly/3 







Lya 
12.132A 


Ly/3 
10.239A 


r 

13.447A 


Lya 
18.967A 


Ly/3 

16.006A 


r 

21.602A 


16.071A 


Ne X 


O VIII 


II Peg 


HEG 


1364±70 


160±28 












8.5 ± 1.6 






MEG 


1217±28 


178±8 


365±21 


2050±70 


419±26 


250±50 


77±16 


6.8 ±0.3 


5.6 ±0.5 


IM Peg 


HEG 


410±30 


80±20 












5.1 ±0.9 






MEG 


385±13 


71±2 


77±11 


530±40 


145±22 


52±25 


55±13 


5.4 ±0.4 


5.0 ± 1.2 


AR Lac 


HEG 


727±30 


106±11 












6.8 ±0.8 






MEG 


632±16 


94±5 


107±14 


820±40 


175±15 


156±44 


121±14 


6.7 ±0.4 


9.9 ±1.9 


HR 1099 


HEG 


2180±40 


288±13 












7.6 ±0.4 






MEG 


1705±21 


222±6 


571±16 


2860±60 


549±19 


400±40 


225±14 


7.7 ±0.2 


7.6 ±0.5 
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Table 3. Plasma densities and Temperatures. 









II Peg 


IM Peg 


AR Lac 


HR 1099 


n c a 
n c a 


Mg XI 
VII 


(10 12 cm- 3 ) 
(10 10 cm- 3 ) 


5-6±S:g 

o n+6.8 
O.Z_ 1 4 


q 9+2.5 

o.z_ 1A 


< 1.8 


i - u -0.8 


rp b 
rp b 


Nc 



(MK) 
(MK) 


6.7-7.3 
4.7-5.5 


7.8-8.8 
4.5-7 


8.1-9.6 

4-5 


6.0-6.4 
4.7-5 



a From He-like triplet diagnostics (Paper I). 
b From Lya/r ratio diagnostics. 
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Table 4. Path length derived from measured Lya/Ly/3. 



Source 


Ion 


Element 3, 


lr 


^RTV 








Abundance 


[cm] 


[cm] 




II Peg 


VIII 


8.97 d 


9.5 • 10 9 


1 • 10 9 


0.04 


IM Peg 


VIII 


9.37° 


1.7- 10 10 


2.2 • 10 9 


0.019 




Ne X [HEG] 


8.86° 


1.6- 10 8 


2.8- 10 7 


0.0002 




Nc X [MEG] 




2.2 • 10 8 


2.8 • 10 7 


0.00018 



a Expressed on the usual spectroscopic logarithmic scale where 
X/H=log(n(X)/n(H))+12, and n(X) is the number density of the element 
X. 

b Loop length from RTV scaling laws: i RT v ~ T 2 /[(IA ■ 10 3 ) 3 • 2n c k B ] 
c Path length as fraction of the stellar radius 
d From Huenemoerder et al. (2001) 

°Scaled from Huenemoerder et al. (2001) II Peg values by the ratio 
of II Peg to IM Peg photospheric metallicities, [Fe/H]n = —0.4 and 
[Fe/H]iM = 0.0, derived by Berdyugina et al. (1998,1999, respectively) 
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5 6 7 8 9 10 11 
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Fig. I. — O VIII (top) and Ne X (bottom) Lya/Ly/3 ratios for the analyzed sources vs. the 
temperature derived from the Lya/He-like r diagnostics. The solid curve is the theoretical 
ratio from the APED database. For the Ne X we show the results from both the HEG 
(filled symbols) and the MEG (empty symbols) measurements; the HEG ratios are shifted 
by +5% on the T axis in order to distinguish the error bars corresponding to HEG and MEG 
measurements for the same source. 
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Fig. 2. — Comparison of the spectra for sources showing optical depth effects and sources 
appearing effectively optically thin. Top - O Vlll Lya, (3 spectral region for II Peg (filled 
profile) and HR 1099 (shifted in A by +0.06A for better readability). Spectra are nor- 
malised such that they have the same unblended Ly/5 line strength, i.e. after correction for 
the Fe XVIII line flux (see text). The different strengths of the Fe line complex near 17 A 
reflects primarily a different Fe/O abundance ratio in the two stars (Drake et al. 2001, Huen- 
emoerder et al. 2001). Bottom - Ne X Lya, (3 spectral region for IM Peg (filled spectrum) 
and AR Lac (shifted in A by +0.03A for better readability). Spectra are again normalized 
to the observed intensities of the Ly/3 line. 
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